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A THEORETICAL ANALYSIS OF THE FREE VIBRATIONS OF 


RING- AND/OR STRINGER-STIFFENED ELLIPTICAL 
CYLINDERS WITH ARBITRARY END CONDITIONS, 

VOLUME I - ANALYTICAL DERIVATION AND APPLICATIONS 

N 

By Donald E. Boyd and C. K. P. Rao 
SUMMARY 

An analysis was made in this study to determine the natural 
frequencies and mode shapes of ring- and/or stringer-stiffened noncir- 
cular cylinders with, arbitrary end conditions. Cases of circular, 
noncircular, unstiffened, and stiffened cylindrical shells with various 
end conditions were investigated and the following observations were 
made. 

1) Comparisons with previous results from experimental and analytical 
studies of circular, noncircular, unstiffened, and stiffened cylindrical 
shells with arbitrary end conditions showed this method of analysis to 
be accurate. 

2) The natural frequencies obtained in this study for a clamped-free 
circular cylinder were slightly higher (for the whole range of m and n) 
than those previously obtained experimentally. This discrepancy 
increases as the number of circumferential waves decreases. 

3) Comparisons with analytical results obtained previously for 
stringer-stiffened, freely supported, circular shells showed that the 
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frequencies previously obtained (neglecting insurface inertias and 
employing Donnell's shell theory) were slightly higher than those of the 
present analysis. The discrepancies between the results decreased as 
the number of circumferential waves increased, which is a typical 
characteristic of Donnell's theory. 

4) Comparisons with Forsberg's exact results for ring-stiffened 
circular shells, showed that the results of the present analysis were in 
error only by a maximum of 0.517> for zero-eccentricity rings and 1.754 
for negative-eccentricity rings. 

5) Comparisons with Al-Najafi and Warbur ton's finite element and 
experimental results (obtained for ring-stiffened circular shells) 
showed that the results for freely supported cylinders obtained from 
the present analysis were closer to their experimental results than 
their analytical results using the finite-element method. For the free- 
free case, of the six experimental results presented, the results of the 
present analysis were closer to the first three experimentally obtained 
frequencies, whereas their finite element results were closer to the 
next three frequencies. 

6) The number of terms required in the displacement series for con- 
vergence of results for ring-stiffened shells differed from problem to 
problem. Shells with positive eccentricities needed more terms for 
convergence than those with zero or negative eccentricities. 
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INTRODUCTION 


Discussion 

The free vibrations of ring- and/or stringer-stiffened circular and 
noncircular cylindrical shells are of interest to designers of flight 
and marine structures. Frequently, fuselages of flight structures and 
hulls of submarines have noncircular cross-section due either to special 
internal storage requirements or to imperfections occurring during man- 
ufacture. The method of analysis developed in this report is capable of 
evaluating the free-vibrational characteristics of ring- and stringer- 
stiffened "singly" symmetric noncircular cylinders with arbitrary end 
conditions . 

Background 

Solutions for the vibrational characteristics of the special cases 
of unstiffened, circular, isotropic cylinders with specialized boundary 
conditions have been available for many years. Recent investigations 
have taken advantage of computers to analyze more complicated models of 
shell structures. One of the most general cases that can be analyzed is 
a stiffened, noncircular, anisotropic cylinder with arbitrary end 
conditions . 

Great attention has been paid to the application of the finite 
element and finite difference methods of analysis because of their 
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generality and adaptability to the computer. However, computer storage 
and the speed of execution are two important factors which have still 
prevented economically feasible studies of shell structures. The closely 
related and well-known Rayleigh-Ritz method was successfully employed in 
the present study to obtain the vibrational characteristics of stiffened, 
noncircular cylinders with arbitrary end conditions. This method may 
provide significant economical advantages over the finite element and 
finite difference methods. The limitation of the Rayleigh-Ritz method is 
that the accuracy of the results is dependent upon the assumed mode 
shapes. In cases such as stiffened, noncircular cylinders with arbitrary 
end conditions (for which the displacement functions can be approximated 
fairly accurately by a double finite series) the Rayleigh-Ritz method is 
certainly advantageous to use. 

Studies of noncircular cylinders are relatively few compared to 
those of circular cylinders. The variable radius of curvature of the 
cylinder causes difficulties in obtaining analytical solutions. If 
finite trigonometric series are used to represent the components of the 
assumed displacement functions, there will be coupling of the circum- 
ferential terms due to noncircularity of the cross-section of the shell. 
Furthermore, the resulting set of simultaneous equations is sufficiently 
large that a digital computer must be used for the solution of the 
general problem. 

Kempner (1) presented energy expressions and differential equations 
for cylindrical shells with arbi trary cross-sections. Ke mpner _and_his_^, 
associates have used these equations to study a wide range of problems 
dealing with statics, buckling and postbuckling (2-7) of a special class 
of oval cylinders. Klosner and Pohle (8, 9, 10) studied the free and 


4 



forced vibrations of the same class of oval cylinders, but with infinite 
length. An approximate method was used in which the frequencies of non- 
circular cylinders were determined by perturbation of the equivalent 
circular cylinder frequencies. Culberson and Boyd (11) obtained exact 
free vibrational characteristics of the same class of oval cylinders 
studied by Klosner and Pohle and showed that the approximate pertur- 
bation technique is accurate for small eccentricities. 

The displacement functions used by Boyd (12) in a static analysis 
of noncircular panels subjected to uniform normal pressures were used in 
a free vibrational analysis of noncircular cylindrical panels by Kurt 
and Boyd (13). 

Herrmann and Mirsky (14) investigated the longitudinal, torsional, 
and flexural vibrations of elliptical cylinders. Malkina (15) also 
studied the free vibrations of oval cylinders. 

Sewall et al. (16, 17) carried out both analytical (by Rayleigh- 
Ritz) and experimental analyses of elliptical unstiffened cylinders with 
arbitrary end conditions. 

Analyses of stiffened shell structures may be classified either as 
"smeared," or as "discrete" depending upon the treatment of the 
stiffeners. In the conventional smearing technique (which is reasonably 
effective if the stiffeners are closely spaced) the effects of the 
stiffeners are averaged out over the entire surface of the shell, thus 
effectively replacing a stiffened shell by an equivalent orthotropic 
shell. A discrete analysis (which is accurate irrespective of the 
number and location of the stiffeners) treats the stiffeners as discrete 
elastic structural elements. 
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The present analysis may be considered as an extension (to include 
noncircularity) of the work in Refs. (l8 and (19) in which the free vibra- 
ational characteristics of ring- and stringer-stiffened noncircular cyl- 
inders with arbitrary end conditions were developed through the use of 
a Rayleigh-Ritz technique. The stiffeners may be arbitrarily located 
and all stiffeners need not possess the same geometric and material 
properties; however, the stiffeners are assumed to be uniform along 
their axes. The analysis considers the extension and flexure of the 
shell and extension, torsion, and flexure about both cross-section axes 
of the stiffeners. The stringers may have nonsymmetric cross-sections 

but the rings are assumed to have "singly" symmetric cross-sections. 

» 

The rotary inertia of the shell is neglected. 

The derivation of the energy expressions for noncircular cylinders 
is described in the Method of Analysis section of this report. The 
stiffener energies are presented in Appendix B. The compatibility 
relations used in these equations are derived in Appendix A. The 
elements of the mass and stiffness matrices are given in Appendix C. 
Documentation of the computer program developed for this analy sis is 
given in Reference (20). 
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SYMBOLS 


A ,A 

H rk 


e x ,e 0 ,e x0 


< e x ) * (e e> 

s r 


E sl’ E rk 
< GJ W < GJ >rk 


IRl^ 

to 

00 
r— 1 

Pi 

M 

IR2J 

to 

IR2 io 

ir3 x 

to 

IR3 2 

IR4 X 

to 

IR4 5 

IR5 1 

to 

IR5 18 

IR6J 

to 

IR6 11 

IX L to IX 5 

ISlj 

to 

IS1 9 

IS2 1 

to 

IS2 

5 

1 yysi 

> I 

xxrk 


yzsi 


length of the shell (except as noted in Figure 1) 
cross-sectional area of the l stringer, k ring 
isotropic plate flexural stiffness 
strains of shell (see eq. (1)) 

normal strains of stringer and ring, respectively 

Young's modulus of shell 

tli thi 

Young's modulus of l ■ stringer, k ring 

til til 

the torsional stiffness of the i stringer, k ring 
thickness of shell 


circumferential integrals of ring equations (see eq. C7) 


longitudinal integrals (see eq. C3) 

circumferential integrals of shell equations 
(see eq. C2) 

th th 

the moment of inertia of the l stringer, k ring 
cross-sectional area, about y' and x axes passing 
through their shear centers 

product of inertia of the stringer cross-sectional 
area about y' and z‘ axes passing through its shear 
center 
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*zzsj l’ ^zzrk 
^yycsA ’^xxcrk 


I 


yzcsX, 


I 


zzcs i 


titx til 

the moment of inertia of the l stringer, k ring 
cross-sectional area about z / axis 

the moment of inertia of the stringer, k^ ring 
cross-sectional area about axes parallel to y and x 
axes passing through its centroid 

til 

product of inertia of the l stringer cross-sectional 
area about axes parallel to y and z axes passing 
through its centroid 

the moment of inertia of the stringer cross- 
sectional area about an axis parallel to z axis 
passing through its centroid 


K 

L 

M* 

N* 


q (t) 

mn ' 
R 


crk 


total number of rings 
total number of stringers 

final value of m in the assumed displacement series 

final value of n in the assumed displacement series 

generalized coordinate 

radius of curvature of the shell 

til 

radius of the centroid of the k ring 
time 


T kinetic energy 

u,v,w longitudinal, circumferential, and radial displacements 

of the middle surface of the shell, respectively 
(see fig. 1) 


u^jV^jW^ displacements of an arbitrary point in the cross- 

1 section of the i 4 -* 1 stiffener in the x, 0, and z 

directions 

u . ,v . ,w_ . displacements of the shear center of the i^ stiffener 
SCI sex SCX , i a j • , . 

in the x, 6, z axrectxons 

u ,v ,w generalized coordinates for symmetric mode displace- 

mn mn mn . . , 

ments u,. v, and w, respectively 

u 7 ,v 7 ,*/ generalized coordinates for antisymmetric mode 

mn’ mn’ mn , _ , ■ , 

displacements u, v, and w, respectively 

U strain energy 
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“*<*> 

W (x) 
m 

x, 0 , z 


axial mode functions representing displacements in the 
x, 8, and z directions 


longitudinal, circumferential, and radial shell 
coordinates (see fig. 1) 


x' ,y' ,z' 


longitudinal, circumferential, and radial coordinates 
of the stiffener, measured from its shear center 


X m< x > 


^lsf, 


"2s Z 


z lsr Z lrk 


Z 2s£’ Z 2rk 


$ (x) 
m 

A -- ,B 
mn,mn mn,mn 

C — ,D 
mn,mn mn ,mn 

E -- ,F 
mn,mn mn ,mn 

N __ ,NN 
mn,mn mn,fflfi 


Bernoulli-Euler beam eigenfunctions 


, th 


y-distance of the shear center of the Z stringer from 
the z axis passing through its point of attachment 

th 

y-distance of the centroid of the stringer from its 
shear center 

th th 

z-distance of the shear center of the Z stringer, k 
ring from the middle surface of the shell 

til tlx 

z-distance of the centroid of the j l stringer, k 

ring from its shear center 

general axial mode function (see eq. 28) 


elements of the stiffness matrix (see appendix C) 


p --,0 -- > elements of the mass matrix (see appendix C) 

mn,mn Tnn,mn i 


R -- , S 
mn,mn mn ,mn 

V 

Pc 

p sjT p rk 

u) 


Poisson's ratio 


mass density of the shell 

th th 

mass density of Z stringer, kiring 
circular frequency 


Subscripts: 

a 

c 


antisymmetric 

refers to cylinder; centroid 
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k 


l 

m, m 

n, n 
r 

s 

sc 

Notes : 


th 

refers to the k ring 
refers to the stringer 

identifies m and m longitudinal modal components 
til til 

identifies n and n circumferential modal components 
refers to rings 
refers to stringers 
refers to shear center 


(1) A comma before a subscript denotes partial differentiation 
with respect to that subscript; 

3 u ^2 w 

e.g., u, x denotes ^ and w >qq denotes 

(2) Superscript T denotes transpose of a matrix. 

(3) Dots over quantities denote differentiation with respect to 
time. 
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METHOD OF ANALYSIS 


The analytical method employed in this analysis was the well-known 
Rayleigh-Ritz (i.e. "assumed modes") energy technique. At the outset 
the strain and kinetic energies of the shell, ring, and stringer were 
derived. The compatibility relations were developed to express the 
displacements of rings and stringers in terms t>f the displacements of 
the median surface of the shell. The total strain energy of the shell 
and that of rings and stringers were combined to obtain the total strain 
energy of the stiffened cylinder expressed in terms of displacements of 
the shell median surface. The total kinetic energy of the stiffened 
cylinder was similarly formulated. Finite series were assumed repre- 
senting the circumferential, -axial, and radial displacements of the 
median surface of the shell and satisfying the shell kinematic boundary 
conditions. Simple trigonometric functions were used to describe the 
circumferential displacement distributions and beam functions were chosen 
to describe distributions along the axis of the shell. The assumed 
displacement functions with undetermined coefficients were substituted 
into the total energy expressions of the structure, and the regular 
eigenvalue problem was formulated by minimizing the action integral. 
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Geometry 


Strain-displacement relations : The classical theories of thin 

shells and beams were used to derive the energy expressions for the 
shell and the stiffeners, respectively. The geometry of the middle sur- 
face of a typical elliptical shell is illustrated by Figure 1. The 
three orthogonal coordinates x, 8, and z locate points within the struc- 
ture and u, v, and w are the corresponding displacement components. The 
variable radius of curvature of the shell cross-section is expressed as 
a function of the 9 coordinate. The following Kempner (1) relations were 
used to determine strains at points within the shell: 


e = u, - zw, 
x x xx 




e = li + j[M£l v ■ z (2R+z) w 

x0 R+z R ’x R(R+z) s x0 

where e , and e Q are normal strains of x- and 0-oriented line elements, 
x y 

respectively, and e is the distortion angle between these two line 

xU 

elements. Furthermore, u, v, w, and R refer to middle surface (z=0) 
values . 

For the stringers and rings the normal strains were expressed as 


(e ) = u 

x s ,x 

s * 


(e e ) r = r 7 ( V r,e + M r) ^ F;( V r ,0 + w r ) 


( 2 ) 


(3) 


'where - the - subscripts - s~ and r _ indicate“arbitrary points-in—the-stringer- 

and ring, respectively. (e ) is the normal strain of the stringer in 

x s 

the x direction, and (e ) is the normal strain of ring in the 0 

0 r 

direction. R is the radius of the centroid of the ring, 
cr 
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Figure 1. Geometry of an Elliptical Shell. 
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Compatibility relations: The geometric details of eccentric 

stiffeners are shown by Figures 2 and 3. The compatibility equations 
relating the displacements of any point in the stiffener cross-section 
to those of its shear center are presented in Appendix A. The following 
equations were derived to determine the displacements in the stiff- 
eners ; - 

For the stringers: u = u - z 7 w - y 7 v (4) 

s scs . scs,x scs,x 


For the rings : 


x z 

v = v - u - — 

r scr R scr,0 R 
scr 


z ( \ 

n - v 

\ scr.Q scr/ 
scr * D 


w = w + x w 
r scr scr,x 


(5) 


where the subscript sc identifies the shear center, and the coordinates 
x 7 , y 7 , and z 7 are measured from the shear center of the 6tiffener. 

The following compatibility equations relating the displacements of 
the shear center of the stiffener to those of the shell's median surface 
were derived and are presented in Appendix A. 

For the stringers: 

u = u - z, w« - y, v, 
scs Is x J ls x 


/ w ’ 


V - V 

scs 


Z ls V 


li 

R " R / 


w scs = w + yis (~r - 1) 


For the rings: 


u = u - z, w, 
scr lr ’x 


/ . Z lr \ Z lr 

V, - l 1 + — r - — 


w = w 
scr 


( 6 ) 


(7) 
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Figure 3. Geometric Details of an Eccentric Stringer Stiffener 
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Strain and Kinetic Energies 


Shell energies : From Reference (1), the strain energy in an isotro 

piCj elastic body subjected to small strains e , e Q , and e Q is 

X O Xo 

u ' I ™pr{ e V 2v Ve + ^ e xe ] d < TOl) <8 > 

vol 2(L_V ' 

For a shell of uniform thickness h, the above expression can be written 


as 


h 

a 2tt 2 


U = 
c 


2(l-v 


- s n [ 

} 0 0 -h 

2 


e x + e2 9 +2ve x e 9 


+ IkVi e 2 

+ 2 x9 


] (R + z) 


dz d9 dx 


(9) 


where E c is Young's modulus and v is Poisson's ratio of the shell. 


After substituting the Equation (1) into the Equation (9) and inte- 
grating over the thickness of the shell, we obtain the strain energy of 
the shell in terms of the displacements of its median surface; i.e. 

u c*TtJ I [ Ru? x + ii ? i (ff + T2 ! k 3 ‘) u? e + 2vu 's v >e 

h o o 

+ (i-v) ■,,».«*•<« + S v ?e + ( R + |s )< 

+ r v - 0 u+ (I + T2T 5 ')" 3 ] d9dx + D J I t' 2 "-*"- 


XX 


0 0 


. (1-V) 2v _ 3(l-v) . „ a . ’90 

+ u 'e"’*9' T ’6 'xx' “ht v -x"'x9 + Ru -xx + — 

(""'ee +w> e8") + s(”'xx“-ee + "-08”-xx) 


211^1 a ] 

R ’x9 J 


d9 dx 
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+ D 




2 

R 



vw, 


e 



+ ”’ee w, e) + jT (D, 0 ("" ; 0 + *-o”) + v (D 

+ •■e’O] 40 dx <10) 

where 

Eh 3 

c 

D = 12(l-v* ) 

The last integral in Equation (10) vanishes for constant R. The first 
two integrals are equivalent to those developed by Miller (21) and by 
Egle and Soder (19). 

Neglecting the contribution of rotary inertia, the shell kinetic 
energy may be written as 

a tt 

T c = p c h J J [u 2 + v 2 + w 2 ] R de dx (11) 

0 0 


where p c is the mass density of the shell and the dot represents the 
time derivative. 

Ring energies : The ring is assumed to be subjected to normal 

strains and shearing strains due to twisting. The cross-section of the 
ring is assumed to be symmetric with respect to the outward normal to 
the shell surface through the line of attachment. The total strain 
e nergy in K rings due to normal strains is 

U r -I ¥j J [< e e> ] 3 “rk « < l2 > 

k-1 0 4,. ' 

rk 
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Using the strain-displacement relation of the ring (Equation (3)) the 
above expression may be written as 

K E 2tt 

“r-I-rJ J F^[''rV"r + V r,e"r + "r V r ,e] “rk d9 < U > 

k=l 0 A . cr X=X k 

rk 

Substituting the first set of compatibility relations of the ring 
(Equations (5)) into Equation (13) and performing the integration over 
the cross-section of the ring, the strain energy of the ring due to 
extension (normal strain) may be written in terms of the displacements 
of its shear center as 

U = U (u , v , w ) (14) 

r r scr scr scr 

ext ext 

The function U (u , v , w ) is given in Appendix B. Com- 

r scr scr scr 

ext 

bining Equations (7) and (14) results in 

U r = U r ^ U ’ V ’ ^ < 15 > 

ext ext 

The function U r t (u,v,w) is also given in Appendix B. 
ext 

The strain energy due to twisting of the rings may be written as 
(Reference 27) 


U r = I 


2rr u 


scr ,9 


w, x6 n 3 


cr x=> 


R d0 
cr 


where (GJ)^ is the torsional stiffness of the k ring. Substitution 

of Equations (7) into Equation (16) results in 

U r = U r (u,v,w) (17) 

tor tor 


The function U r (u,v,w) is given in Appendix B. 
tor 

The kinetic energy of the ring is 
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T r - i 1 P rk J I [ S r + \ + K ] “rk *c 
k=l 0 A - 


Substitution of Equations (5) into the above equation and integrating 

over the cross-section of the rings, and then substituting the Equations 

(7) into the resulting expression we have, 

T = T (u, v, w) (19) 

r r 

The function T^Cu, v, w) is given in Appendix B. Note that Equation (19) 
includes both translation and rotation effects. 

Stringer energies : The stringer is assumed to be subjected to both 

extension and twisting. The cross-section of the stringer may be non- 
symmetric. The strain energy due to normal strain in the stringer is 


V -X-rJ I [ ( \>J “ 


4 = 1 0 A 


8 " 0 = 0 , 


or, introducing Equation (2), 


u =f fair j r u r dA dx 

s L 2 J J L s,x J s£ 


1=1 0 A 


Substitution of Equation (4) into the above equation and integrating 

over the cross-section of the stringer, and then substituting 

Equations (6) into the resulting expression we obtain 

U = U (u,v,w) (22) 

ext ext 

The function U s (u,v,w) is given in Appendix B. 

_ — ext— — — ■ — ; — ^ 

The strain energy due to twisting of the stringer may be written as 
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where (GJ) . is the torsional stiffness of the stringer. Thus, 
sZ 


L a 2 2 

(GJ) sA „ w, „ v. 


* -7 ^stf + 

s L> 2 J L t,3 _2 t>2 J 


dx 


tor A=1 0 R R 

The kinetic energy of stringer is 
L a 


J e=e ( 


sJL 

combining Equations (4, 6, and 25) and integrating the resulting 

expression over the cross-section of the stringer results in 

T = T (u,v,w) 
s s 


The function T g (u, v, w) is given in Appendix B. 


( 24 ) 


f S J I [ a s 3 + 1 s 2 +i sl. “.J dx (25 > 

i-1 0 A , e=0 A 


(26) 


Displacement Functions , 

The displacements, u, v, and w were assumed to be double finite 
series. Each term of the series is a product of a circumferential and 
an axial modal function weighted by a time-dependent generalized 
coordinate (unknown amplitude coefficient). The assumed displacement 
functions were: 

M* N* 

u(x, 9, t) = / / (u cos n0 + u' sin n0) U (x) e' UOt 

' Aj mn mn m v 

m=0 n=0 


M* N* 


v(x, 0, t) I (v^ sin ne - v'^ cos n0) Vjx) e 


hot 


m=0 n=0 


M* N* 


w(x, 0, t) =V Y (w cos n0 + w ' sin n0) W (x) e^*" (27) 

Li Lt mn mn m 


m=0 n=0 
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where U (x) , V (x) , and W (x') are the axial mode functions which satisfy 
m ni in J 

at least the kinematic boundary conditions of the stiffened shell. Also, 

u , v , and w are unknown amplitude coefficients of the symmetric 
mn mn mn 

circumferential modes, and u 7 , v' , and w* are those associated 

* mn mn mn 

with the antisymmetric modes. 

In this analysis the axial mode functions U (x) , V (x) , W (x) were 

m m m 

expressed by a single function 4^(x) such that 

u (X) = ij- 4 (x) (28a) 

m dx m 

V(x) = # (x) 

m m 

W (x) = 4 (x) 

m m 

The following functions were implemented in this analysis. 


Boundary Condition 

Function Used 

Eqn. No. 

Freely supported: 

§ (x) = \IT Sin — 
m a 

(28b) 

Clamped- free: 

4 (x) = X , Characteristic function of 
m m 

a Clamped-free beam. 

(28c) 

Clamped-clamped: 

4 (x) = X , Characteristic function of 
m m 

a Clamped-clamped beam. 

(28d) 

Free-free: 

§ Q (x) = 1 

V*)=!-s 

4 (x) = X . , Characteristic function 
m 7 m-1 

(28e) 


of a Free-free beam, (ms 2) 



The characteristic functions X , their derivatives and eigenvalue 

m 

properties are tabulated in Reference (22) . 
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The Frequency Equation 


The total strain energy of the stiffened shell was obtained by 


combining Equations (10, 15, 17, 22, and 24). Similarly, the total 


kinetic energy was obtained by combining Equations (11, 19, and 26). 


Substituting Equations (27 and 28) into the total energies of the 


stiffened shell, the strain energy expression becomes a positive definite 


quadratic function of the generalized coordinates u , v , w , u' , 

ton mn nm mn 

v' , and w* . Furthermore the kinetic energy expression becomes a 
mn mn 

positive definite quadratic function of the generalized velocities u , 
‘ mn 

O o O / O / O / 

V , W , U , V , w 
mn mn mn mn mn 

The total strain energy of the structure may be written as 


M* N* M* N* 


total 


I! E ^mn,mn Snn'^mn 


m=0 n=0 m=0 n=0 


(29) 


where 


d 2 U 


total 


total 


S q Bq-- dg Bcr~ 

inn mn mn mn 


= K __ = K„_ 
mn ,mn mn,mn 


are known as elements of the stiffness matrix. 

The total kinetic energy of the structure may be written as 

M* N* M* N* 

T total ■*I l l mn ,inn Srm^mn 

m=0 n=0 ffl=0 n=0 


(30) 


where M __ are the elements of the mass matrix, 
mn ,mn 

The mass and stiffness matrices obtained by the above operations 
were used together with Hamilton's principle to formulate the regular 
eigenvalue problem resulting in 
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q S ) 

= 0 (31) 

q a ) 

where K, and M represent stiffness and mass matrices of size 
3(M*+1) (N*+l) , q and q denote the symmetric and antisymmetric mode 

S 3. 

vectors, respectively, and superscript T denotes the transpose of a 
matrix. 

In Equation (31) the off-diagonal submatrices of both the 
stiffness and mass matrices vanish if the cross-section of the 
stiffened shell is symmetric with respect to the vertical axis (where 
9=0). Thus, the above equation is uncoupled into two equations; one 
for symmetric, and the other for antisymmetric modes. The equation for 
the symmetric mode problem may be written as. 


’ K 

ss 

K 

sa 

3 

' M 

ss 

M 

sa 

T 

K 

K 

" U) 

T 

M 

M 

. sa 

aa _ 


sa 

aa _ 



Each letter in the stiffness and mass matrices represents a submatrix 
(presented in Appendix C) of order (M*+l) (N*+l) . 
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COMPUTER SOLUTION 


General 

A computer program (20) was developed to find the eigenvalues and 
eigenvectors of Equation (32). The mass and stiffness matrices were 
generated in this program and the frequencies and mode shapes were 
computed using the subroutine EIGENP (23). The Oklahoma State University 
IBM Model 360/65 computer was employed for this project. 

The input data to the program may be categorized into four kinds. 

The first kind is general data. For example, the title of the problem, 
number of terms considered in the assumed displacement series, whether 
or not the cross-section of the shell is circular, the number of 
stiffeners, etc. The other three kinds of data are shell data, stringer 
data, and ring data. 

The radius of curvature (R) of the shell was considered to be a 

function of the 0-coordinate. The expressions for R, (■£■),, and (R) , , 

K 0 0 

were calculated (considering elliptical cross-section) in the function 
subprograms (RSHL) , (RRRT) , and (RSHLT) , respectively. This procedure 
was used to make the computer program capable of analyzing arbitrary 
singly symmetric stiffened oval cylinders. However, only elliptical 
cylinders were considered in the present study. 
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Natural Frequencies and Mode Shapes 


If the number of circumferential and axial terms considered in the 
assumed displacement series are M and N , respectively, (including 
m = 0, and n = 0, when needed) then the order of the stiffness and mass 
matrices is 3M N . Equation (30) may be written as 



K= Stiffness Matrix / 

M = Mass Matrix 

ou = The natural frequencies from : Equation (33) in radian/sec. 


If the matrices K and M became singular due to the presence of 
zeroes in some of the rows and columns, the matrices were condensed by 
eliminating those rows and columns of zeroes. The subroutine called 
EIGENP (23), with double precision, was used to calculate the frequencies 
(uj ^ of Equation (33) and the resulting eigenvectors 
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w 


Once the eigenvalues and eigenvectors were obtained, the corresponding 
mode shapes were found. 
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NUMERICAL RESULTS 


Introduction 

The analysis described in this report was substantiated by com- 
paring the results of this study with some of those obtained by previous 
investigators. Some parametric studies of stiffened noncircular 
cylinders were made and are also presented in this chapter. 

Comparison With Known Solutions 

This section presents the comparison of natural frequencies for 
(1) an unstiffened circular cylinder with various boundary conditions; (2) 
ring- and/or stringer-stiffened circular cylinders with various end 
conditions; (3) unstiffened noncircular shells with various end 
conditions; and, (4) ring- and stringer-stiffened elliptical cylinders. 

Comparison of results for the unstiffened circular shells : 

Forsberg (24) presented exact frequencies for a freely supported 
unstiffened circular cylinder, obtained by solving the differential 
equations of motion. Hie results of this analysis and those of 
Forsberg' s exact solution are compared in Table I. Both the analyses 

used - the— Fltigge - shell— theory.— As— is-evident_from__the_Table_JI,_.gqqd 

correlation exists between the frequencies of both the analyses. Such 
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TABLE I 

COMPARISON OF ANALYTICAL FREQUENCIES OF A FREELY 
SUPPORTED UNSTIFFENED CIRCULAR CYLINDER, 
OBTAINED BY THE PRESENT ANALYSIS 
AND FORSBERG (Hz.) 


n 

m 

PRESENT 

ANALYSIS 

FORSBERG a 


i 

778 

778 

2 

2 

2449 

2449 


3 

4253 

4253 


1 

628 

627 

3 

2 

1458 

1458 


3 

2682 

2681 


I 

974 

974 


2 

1304 

1303 

4 

3 

2021 

2020 


4 

2947 

2946 


a) Reference (24), figure 3(a). 
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type of accuracy was expected because the assumed mode functions 
satisfy the freely supported boundary condition exactly. 

Comparisons were also made with the results of Reference (16) for 
the same boundary condition and m = 1 and 2. These are presented in 
Table II. In Reference (16), Sewall et al., using Sander's shell theory 
(25), applied the Rayleigh-Ritz method as in our analysis. As is evident 
from Table II, excellent comparisons were obtained. 

Figure 4 shows a comparison between the analytical and experimental 
results of Reference (17) and those of the present analysis (for m = 1) 
considering a clamped-free, unstiffened, circular shell. The frequency 
curves reveal that this analysis yields results similar to those of 
Reference (17). The slight differences might be attributed to the 
difference in the shell theories. Comparisons were also made with the 
experimental results of Park, A. C. et al., (26) and the analytical 
results of Egle and Soder (19). These are presented in Table III. In 
this comparison four-place accuracy was obtained between the analytical 
results of Egle and Soder and the present analysis. The discrepancy 
between the analytical and experimental results increases as the number 
of circumferential waves decreases. Egle and Soder speculated in Refer- 
ence (19) that the shell end may not have been absolutely fixed in the 
experiments. 

The experimental and analytical results of Reference (16) for free- 
free circular shells were used to establish the validity of the present 
analysis for this boundary-condition case. Table IV shows the com- 
parison of the results for-m = 1 and-2. - The present analysis yielded 
four-place accuracy. 
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TABLE II 


COMPARISON OF ANALYTICAL FREQUENCIES OF A FREELY 
SUPPORTED UNSTIFFENED CIRCULAR CYLINDER 3 , 
OBTAINED BY THE PRESENT ANALYSIS 
AND SEWALL (Hz.) 



m = 

1 

m 

= 2 

n 

PRESENT 

ANALYSIS 

SEWALL 
(Ref 16) 

PRESENT 

ANALYSIS 

SEWALL 
(Ref 16) 

i 

1565.3 

1565.0 

2309.3 

2309.0 

2 

894.1 

894.1 

1782.4 

1782.0 

3 

529.8 

529.8 

1314.9 

1315.0 

4 

338.6 

338.6 

968.4 

968.4 

5 

235.6 

235.6 

726.3 

726.3 

6 

182.1 

182.1 

560.3 

560.3 

7 

162.2 

162.2 

448.6 

448.6 

8 

166.9 

166.9 

377.2 

377.2 

9 

188.6 

188.6 

338.1 

338.1 

10 

221.3 

221.3 

325.7 

325.1 

11 

261.7 

261.7 

335.0 

335.0 

12 

308.0 

308.0 

361.0 

361.0 

13 

359.5 

359.5 

399.6 

399.5 

14 

415.6 

415.6 

447.5 

447.5 

a) The 

geometry of 

the shell 

is given in 

Reference (16) 
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TABLE III 


COMPARISON OF ANALYTICAL AND EXPERIMENTAL FREQUENCIES 
OF A CLAMPED-FREE UNSTIFFENED CIRCULAR CYLINDER 3 

(Hz.) 




m = 1 



m = 2 


n 

EGLE & b 

PRESENT b 

PARKS C 

EGLE & b 

PRESENT b 

PARKS C 


SODER 

ANALYSIS 

et al . 

SODER 

ANALYSIS 

et al . 

2 

104.4 

104.4 

87.2. 
95. 1 & 

- 

508.2 

- 

3 

55.6 

55.6 

51.5 

- 

281.3 

- 

4 

52.0 

52.0 

50.4 

177.9 

177.9 

168. 5 X 
170. 2 & 

5 

- 

71.6 

70.9 

- 

135.4 

132.8 

6 

- 

101.8 

101.4 

- 

132.0 

128. 8 r 
130. 1 & 

7 

139.1 

139.1 

138.8 

154.2 

154.2 

153.6 

8 

182.6 

182.6 

182.2 

191.2 

191.2 

191.3 


a) Reference (19), configuration 1, p. 28. 

b) Fliigge shell theory, insurface inertias included. 

c) Reference (26), model 1. 
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TABLE IV 


COMPARISON OF ANALYTICAL AND EXPERIMENTAL 
FREQUENCIES OF A FREE-FREE UNSTIFFENED 
CIRCULAR CYLINDER (Hz.) 




m = 1 



m = 2 


n 

PRESENT 3 

SEWALL b 

SEWALL b 

PRESENT 

SEWALL 

SEWALL 


ANALYSIS 

ANALYSIS 

EXPERIMENT 

ANALYSIS 

ANALYSIS 

iXPERIMENT 

i 

2012. 0 C 

2014. 0= 

- 

2288.0 

2293.0 

- 

2 

7.5 

7.5 

7.7 

1613.0 

1616.0 

- 

3 

19.0 

19.0 

18.9 

1066.0 

1068.0 

- 

4 

34.2 

34.2 

35.7 

716.9 

717.8 

- 

5 

53.4 

53.4 

53.0 

504.4 

504.8 

- 

6 

76.6 

76.7 

76.4 

375.4 

375.6 

377.3 

7 

104.1 

104.1 

103.8 

299.8 

299.9 

299.1 

8 

135.7 

135.7 

135.3 

262.6 

262.2 

257 .4 
262. 1 & 

9 

171.4 

171.5 

170.7 

253.6 

253.4 

248 . 8 
249. 3 & 

10 

211.4 

211.5 

210.2 

266.5 

266.3 

268.8 

11 

255.6 

255.7 

253.0 

294.8 

294.7 

290,9 

12 

303.9 

304.1 

305.5 

333.9 

334.0 

327.6 

13 

356.5 

356.7 

352.0 

381.2 

381.1 

- 

14 

413.3 

413.5 

412.5 

434.7 

434.7 

436.6 


a) Fliigge shell theory; 6 even, and 6 odd axial mode functions 
considered. 

b) Reference ( 16) . 

c) Extensions 1 frequency. 
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Comparison of results for stringer-stiffened circular shells : 

Egle and Sewall (18) presented frequencies obtained for a stringer- 
stiffened, freely supported, circular cylinder using a method similar to 
that of the present analysis but using the Donnell shell theory and 
neglecting the insurface inertias of the stiffened shell. Hie shell 
theory used in the present analysis was modified to Donnell theory in 
order to compare the results of this analysis with those of Egle and 
Sewall. Table V gives the comparison between the frequencies for m = 2. 
The frequencies of Egle and Sewall are slightly higher than those of the 
present analysis, evidently attributable to their neglect of the inplane 
inertias. It is evident from Table V that the discrepancy between the 
results of both the theories decreases as the number of circumferential 
waves increases, which is a typical characteristic of Donnell theory. 

Comparison of results with ring-stiffened circular shells : 

Forsberg (24) obtained exact solutions for the natural frequencies of 
ring-stiffened circular cylinders. Bushnell (27) obtained the natural 
frequencies of ring-stiffened segmented shells of revolution using an 
energy method in conjunction with the method of finite differences. 

The compatibility relations and the energy expressions used by Bushnell 
are similar to those of the present analysis. Table VI presents the 
frequencies obtained by Forsberg, Bushnell, and the present analysis for 
freely supported circular cylinders with three rings of both zero and 
negative eccentricity. The frequencies of this analysis which are 
presented in Table VI were obtained by considering 12 even and 13 odd 
axial mode functions in the assumed displacement series. The results of 
this analysis are in excellent agreement with the exact frequencies 
obtained by Forsberg and the approximate frequencies of Bushnell. The 
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COMPARISONS OF FREQUENCIES OF FREELY SUPPORTED CYLINDERS 
WITH AND WITHOUT INSURFACE INERTIAS 
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TABLE VI 


COMPARISON OP FREQUENCIES OF A FREELY SUPPORTED CIRCULAR 
CYLINDER 3 WITH THREE SYMMETRIC AND INTERNAL RING 
STIFFENERS, OBTAINED BY THE PRESENT ANALYSIS, 
BUSHNELL, AND FORSBERG (Hz.) 





SYMMETRIC 



INTERNAL 


n 

ID 

FORSBERG b 

BUSHNELL C 

PRESENT^ 

ANALYSIS 

FORSBERG 

BUSHNELL 

PRESENT 

ANALYSIS 


l 

788 

787 

787 

999 

987 

994 

2 

2 

2219 

2219 

2219 . 

2254 

2264 

2252 


3 

3796 

3802 

3801 

3710 

3741 

3711 


1 

1155 

1152 

1152 

2087 

2066 

2081 

3 

2 

1661 

1660 

1660 

2397 

2382 

2386 


3 

2617 

2619 

2618 

3073 

3068 

3066 


1 

1988 

1982 

1988 

3161 

3120 

3142 

4 

2 

2132 

2130 

2141 

3085 

3023 

3032 


3 

2535 

2539 

2548 

3014 

3019 

3030 


a) Reference (24), figure 3(a). 

b) Exact solution obtained by solving the equations of equilibrium. 

c) Reference (27), an energy formulation is used in conjunction 
with the method of finite differences. 

d) Energy expressions of ring are similar to those of Reference (27) 
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max imum discrepancy encountered for the case of zero eccentricity ring 
stiffener was 0.51% and 1.75% for the negative eccentricity, ring- 
stiffened case. The external ring-stiffened shell of Forsberg was also 
studied but the frequencies obtained did not converge for 12 even and 
13 odd axial mode functions in the assumed displacement series; hence 
those results are not presented in this report. 

Comparisons were also made with some of the results of Al-Najafi 
and Warburton (28), for freely supported and free-free ring-stiffened 
circular shells and are presented in Table VII. Their results were 
obtained using a finite element technique employing five elements per 
bay. Significant reduction in the order of the matrices was obtained in 
their study by considering the symmetry of the structures and neglecting 
insurface inertias. The results of the present analysis given in Table 
VII were obtained by considering circumferentially symmetric and 10 
even and 10 odd axial mode functions in the assumed displacement series 
but including ihsurf ace inertias. The values for the frequencies con- 
verged for 15 even and 15 odd terms but the difference between the 
results for 10 terms and 15 terms was rather small. Hence, in order to 
compare on the basis of the order of the matrices, the results of 10 
terms was chosen for comparison. It is evident from Table VII that the 
frequencies of the present analysis for the freely supported case are 
lower than those of the finite element method (except for m = 3) and 
are also closer to the experimental values. For the free-free case, 
the finite element results were observed to be closer to experimental 
values than the results of the present analysis, except for m = 1 and 2. 
In general, the agreement between the results of this analysis and those 
of the finite element and the experiment is good. 
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TABLE VII 

COMPARISON OF FREQUENCIES OF RING -STIFFENED CYLINDERS, 
OBTAINED BY RAYLEIGH-RITZ AND FINITE 
ELEMENT METHODS (Hz.) 

(n = 4) ; d = 0.25 in. 


m 

FREELY SUPPORTED 



FREE-FREE 


RAYLEIGH- 3 

RITZ 

FINITE^ 

ELEMENT 

EXPRTL , 

m 

RAYLEIGH- 

RITZ 

FINITE 

ELEMENT 

EXPRTL. 

i 

1867 

1873 

1867 

o c 

1550 

1547 

1551 

2 

2089 

2091 

2076 

l c 

1538 

1537 : 

1539 

3 

2651 

2650 

2600 

2 

1889 

1895 

1890 

4 

3415 

3429 

3355 

3 

2303 

2290 

. 2287 

5 

4239 

4270 

- 

4 

3075 

3044 

3044 

6 

4925 

5022 

- 

5 

3955 

. 3920 

3916 

7 

5846 



6 

4910 

■ - 

- . 

8 

6585 

- 

- 

7 

■ 5548 

- 

- 

9 

7330 


- 

8 

■ 6349 

• • • - : • 

■ - 

10 

8079 

. - 


9 

7103 

- ' 

- 


a) Present Analysis, number of terms considered in the displace- 
ment series is 10. 

b) Reference (28). 

c) Rigid body modes. 






In order to show the rate of convergence of the results of this 
study, the frequencies were obtained with different assumed numbers of 
terms. These results are presented in Tables VIII and IX for the 
freely supported and free-free ring-stiffened shells studied by Al- 
Najafi and Warburton. Tables VIII and IX show that the rate of con- 
vergence of frequencies is rather rapid. 

Comparison of results with ring- and stringer-stiffened circular 
shells : Park, A. G. et al. (26), presented a considerable amount of 

experimental information on the frequencies and mode shapes of stif fened 
and unstiffened circular and elliptical shells with clamped-free ends. 
Egle and Soder (19) compared their analytical results with those of 
Park's experimental results for a clamped-free circular cylinder with 
three equally spaced internal rings and sixteen internal stringers. 

The same shell was analyzed by the present analysis and comparisons are 
indicated in Table X. Because the cross-section of the stiffened shell 
was symmetric with respect to both the vertical and horizontal axes, 
the frequencies of even and odd circumferential modes were able to be 
evaluated separately. It is interesting to notice in Table X that the 
results of the present analysis are consistently lower than those of 
Egle and Soder. This improvement in the frequencies may be attributed 
to the improved stiffener theories of the present analysis. The fact 
that the discrepancy between the analytical and experimental 
frequencies decreases with the increase in wave numbers n and m 
suggests that the boundary conditions of the experiment and the theory 
may not match. 

The results of the present analysis were obtained with 10 axial 
mode functions and 3 even and 3 odd circumferential mode functions. 
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TABLE VIII 


SPEED OF CONVERGENCE OF FREQUENCIES OF FREELY 
SUPPORTED RING-STIFFENED CIRCULAR CYLINDER 3 
(Hz.) , n=4 


v 

c\ 

m \ 

5 

10 

12 

14 

15 

i 

2032.29 

1867.32 

1853.29 

1841.82 

1841.83 

2 

2136.32 

2089.33 

2076.62 

2067.81 

2067.81 

3 

2682.82 

2651.32 

2640.59 

2634.31 

2634.30 

4 

3446.09 

3414.67 

3414.65 

3409.95 

3409.94 

5 

4263.22 

4239.00 

4238.98 

4238.97 

4235.32 

6 

4924.91 

4924.59 

4924 . 58 

4924.57 

4924.47 

7 

5877.52 

5845.98 

5845.97 

5845.97 

5841.54 

8 

6613.81 

6585.41 

6585.39 

6580.90 

6580.89 

9 

7348.25 

7329.87 

7321.17 

7316.42 

7316.41 

10 

8098.23 

8079.40 

8072.23 

8067.25 

8067.24 


a) Reference (28), figure 2(c). 

b) Number of terms considered in the displacement series. 

c) Axial wave number. 
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TABLE IX 


SPEED OF CONVERGENCE OF FREQUENCIES OF FREE-FREE 
RING-STIFFENED CIRCULAR CYLINDER 3 
(Hz. ) , n = 4 


V M b 

c\ 
m \ 

5 

k 

10 

12 

14 

15 

i* 

1591.53 

1549.60 

1546.82 

1546.13 

1544.91 

2* 

1585.73 

1538.16 

1537.45 

1536.33 

1535.35 

3 

2046.65 

1888.92 

1823.09 

1816.19 

1816.05 

4 

2380.46 

2303.22 

2300.84 

2299.44 

2299.35 

5 

3127.52 

3075.50 

3067.22 

3066.92 

3066.66 

6 

3979.47 

3955.27 

3952.06 

3951.22 

3950.53 

7 

4973.26 

4909.71 

4836.28 

4833.91 

4833.57 

8 

5595.02 

5548.42 

5542.69 

5540.21 

5539.64 

9 

6439.71 

6348.83 

6312.89 

6309.63 

6308.67 

10 

7189.93 

7102.58 

7096.81 

7093.99 

7091.25 


* Rigid body modes . 

a) Reference (28), figure 2(c). 

b) Number of terms considered in the displacement series. 

c) Axial wave number. 
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TABLE X 


COMPARISON OF ANALYTICAL AND EXPERIMENTAL 
FREQUENCIES OF A CLAMPED-FREE 
RING- AND STRINGER-STIFFENED 
CIRCULAR CYLINDER (Hz.) 


n 

m 

PRESENT 3 

EGLE & b 

PARK C 

ANALYSIS 

SODER 

et al. 


i 

100.2 

105.8 

80.2 
88. 2 & 

2 

2 

A32.2 

A33.9 

- 


3 

907.0 

- 

- 


1 

207.6 

216.9 

18A.6 

A 

2 

276.0 

285.9 

251.5 


3 

A37.2 

AA7.1 

397. 0 
A30.A & 


1 

308.5 

315.0 

- 

6 

2 

3A5.9 

353.8 

- 


3 

A02.6 

AlA.O 

- 

a) n 

= 2, 

A, 6; m = 1 

to 10. 



b) Reference (19). 

c) Reference (26), model IS. 


A3 





The reason for considering fewer number of circumferential terms than 
the axial terms is that the coupling between the circumferential mode 
functions (due to the presence of stringers) is rather weak. This was 
also noticed experimentally by Scruggs et al. (29). The coupling 
between the axial mode functions (due to the presence of rings) is 
considerable; hence 10 terms were considered in the longitudinal 
direction. To determine whether or not 10 terms were sufficient for 
obtaining reasonably well-converged frequencies, M* was increased to 
30 and only one circumf erential term was used. The comparison between 
these results is shown in Table XI. Since the difference in the results 
was found to be negligible, it was concluded that 10 terms were suffi- 
cient for convergence. 

Comparison of results with unstiffened noncircular shells : Having 

established satisfactory results for stiffened and unstiffened circular 
shells of arbitrary end conditions, comparisons were then made for 
unstiffened noncircular shells. Sewall et al. (16, 17) presented 
analytical and experimental results for elliptical shells with arbitrary 
end conditions. Tables XII and XIII compare the analytical symmetric 
and antisymmetric frequencies for freely supported elliptical shells of 
eccentricities of 0.526 and 0.760 for m = 1. It is evident from Tables 
XII and XIII that the agreement between the results of both Sewall and 
the present analysis is generally satisfactory and is excellent for n 
less than 10. 

Comparison of results obtained for elliptical shells with free- 
free and clamped-free end conditions were also made and are presented 
in Tables XIV and XV, respectively. The results of this analysis 
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TABLE XI 


CONVERGENCE OF FREQUENCIES OF CLAMPED-FREE RING- AND 
STRINGER-STIFFENED CIRCULAR CYLINDER (Hz.) 
(Circumferentially Symmetric) 


n 

m 

a 

b 


i 

99.32 

100.19 

2 

2 

428.66 

432.19 


. 3 

903.77 

906.96 


a) 

k 

N = 2, 

k 

M = 30 

b) 

N* = 6, 

o 

ii 
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TABLE XII 


COMPARISON OF ANALYTICAL FREQUENCIES OF FREELY 
SUPPORTED ELLIPTICAL CYLINDERS 3 (Hz.) 


n 


e = 0.526 

, m = 1 


SYMMETRIC 

ANTISYMMETRIC 

PRESENT b 

ANALYSIS 

SEWALL C 

PRESENT 

ANALYSIS 

SEWALL 

0 

2550.2 

2550.0 

- 

- 

1 

1439.7 

1440.0 

1685.7 

1686.0 

2 

876.6 

876.6 

888.9 

888.9 

3 

524.1 

524.1 

524.2 

524.2 

4 

335.5 

335.5 

335.6 

335.5 

5 

234.3 

234.3 

234.3 

234.2 

6 

184.2 

184.2 

184.2 

184.2 

7 

156 . 9 

157.1 

156.9 

157.0 

8 

160.1 

160.2 

160.2 

160.2 

9 

189.7 

189.8 

189.4 

189.8 

10 

221.5 

221.9 

221.8 

221.9 

11 

260.8 

261.9 

261.7 

261.9 

12 

307.6 

308.1 

307.9 

308.1 

13 

348.9 

359.5 

355.8 

359.5 

14 

405.7 

415.6 

413.9 

415.6 

a) The geometric 

and material 

properties 

of the shells 


are given in Reference (16). 

b) Number of terms used is 13. 

c) Reference (16). 
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TABLE XIII 

COMPARISON OF ANALYTICAL FREQUENCIES OF FREELY 
SUPPORTED ELLIPTICAL CYLINDERS 3 (Hz.) 




e = 0.760 

, m = 1 


n 

SYMMETRIC 

ANTISYMMETRIC 


PRESENT^ 

ANALYSIS 

SEWALL C 

PRESENT 

ANALYSIS 

SEWALL 

0 

2611.8 

2612.0 

- 

- 

1 

1237.7 

1238.0 

1855.7 

1856.0 

2 

785.1 

785.2 

858.5 

858.5 

3 

491.1 

491.1 

492.5 

492.4 

4 

319.8 

319.4 

318.9 

319.4 

5 

- 

- 

226.6 

226.9 

6 

- 

- 

- 

- 

7 

138.5 

138.5 

138.5 

138.5 


140.0 

140.1 

140.1 

140.1 

8 

& 

& 

& 

& 


177.8 

- 178.3 

178.5 

178.3 

9 

182.3 

& 

184.1 

& 

184.0 

184.1 


226.1 

226.9 

- 

- 

10 

221.7 

223.9 

223.5 

223.9 

11 

261.6 

263.6 

259.2 

263.6 

12 

310.6 

307.3 

296.9 

307.3 

13 

378.4 

359.4 

338.6 

359.4 

. 14 

464.8 

, 417.1 

399.6 

417.1 

a) The 

geometric 

and material 

properties 

of the 


shells are given in Reference (16). 

b) Number of terms used is 13. 

c) Reference (16). 
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COMPARISON OF ANALYTICAL AND EXPERIMENTAL FREQUENCIES OF A FREE-FREE ELLIPTICAL CYLINDER (Hz. 
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b) Reference (16). c) Reference (9). d) Reference (10). 



TABLE XV 


COMPARISON OF ANALYTICAL AND EXPERIMENTAL 
FREQUENCIES OF A CLAMPED-FREE 
ELLIPTICAL CYLINDER (Hz.) 




£ 

i = 12.95, b 
m = 1 

= 11.01 





SYMMETRIC 



ANTISYMMETRIC 

n 

PRESENT 3 

SEWALL 

SEWALL 

PRESENT 

SEWALL 

SEWALL 


ANALYSIS 

ANALYSIS 

EXPERIMENT 

ANALYSIS 

ANALYSIS 

EXPERIMENT 

i 

736. 6 

739.2 

- 

838.0 

840.1 

- 

2 

384. 2 

390.6 

- 

387.6 

394.1 

- 

3 

212.4 

217.5 

201.9 

& 

212.4 

217.5 

204.8 


- 

- 

201.1 

- 

- 

- 

4 

133.6 

136.4 

129.5 

& 

133.6 

- 

134.0 


- 

- 

129.1 

- 

- 

- 

5 

97.5 

99.5 

96.4 

97.5 

99.5 

100.2 

6 

94.6 

95.9 

94.2 

& 

94.6 

95.9 

94.5 


- 

- 

93.1 

- 

- 


7 

113.2 

114.2 

115.1 

113.2 

114.2 

116.5 

8 

138.4 

139.6 

141.8 

& 

- 

136.4 

142.3 


- 

- 

140.6 

138.4 

139.6 

- 

9 

171. 2 

171.4 

17 6.0 

171. 2 

171.4 

176.2 

10 

210.0 

210.1 

217.2 

& 

209.9 

210.1 

216.3 


- 

- 

217 .1 

- 

- 

- 

11 

253.6 

253.7 

260.4 

253.5 

253.7 

260.8 

12 

301.4 

301.7 

309.5 

301.6 

301.7 

310.6 

13 

353.8 

354.1 

365.0 

354. 9 

354.1 

• 

14 

409.8 

410.7 

423.6 

412.8 

410.7 

_ 

a) 

N* = 20, M* 

= 2. 
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are similar to those obtained analytically by Sewall. Also included are 
Sewall's experimental results and analytical results obtained by Klosner 
(9, 10). 

Comparison of results with ring- and stringer-stiffened elliptical 
shells: Park, A. C. et al. (26) presented experimental frequencies and 
mode shapes for a clamped-free elliptical cylinder with four equally 
spaced internal rings and sixteen internal stringers. This shell was 

V . , 

also analyzed by the present analysis, and some comparisons are pre- 
sented in Table XVI. Due to the symmetry of the cross-section with 
respect to both the vertical and horizontal axes, the frequencies of 
even and odd circumferential modes were evaluated separately. As is 
evident form Table XVI, the theoretical results are consistently 
higher than the experimental results. The discrepancy between 
the analytical and experimental frequencies may again be attributed to 
the possible difference in the boundary conditions of the experiment 
and the theory. However, storage limitations of the IBM 360/65 computer 
prevented the consideration of a sufficient number of terms in the dis- 
placement series to assure convergence of frequencies. The results of 
the present analysis were obtained with 5 axial mode functions and 6 
even and 6 odd circumferential mode functions. 

Studies of Stiffened Noncircular Cylinders 

Having obtained satisfactory comparisons with known solutions of 
the circular, noncircular, unstiffened, and stiffened cylindrical shells, 
two studies of stiffened noncircular shells were made. This section 
presents the results of those studies. 
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3 


TABLE I 


COMPARISON OF ANALYTICAL AND EXPERIMENTAL FREQUENCIES OF 
A CLAMPED-FREE ELLIPTICAL CYLINDER 3 WITH FOUR RINGS 
AND TWELVE STRINGERS 


n 

m = 

1 

m = 

2 

PRESENT 15 

ANALYSIS 

PARK C 

PRESENT . 
ANALYSIS 

PARK 

i 

177.92 

163.5 

- 

- 

2 

92.08 

60.8, 
79. 7 & 

- • 

- . 

3 

151.75 

141.1 

242.64 

226.7 

4 

- ■ 

- 

377.68 

352.6 


a) The geometry jaf the stiffened shell is 
given in figure 32, model 4S, Park, A. C. 
et al. , dynamics of shell-like lifting 
bodies. Part II, the experimental investi- 
gation. AFFDL-TR-65-17 , Part II, June, 1965. 

b) Rayleigh-Ritz method N* = 12, M* = 5. 

c) Experimental results. 
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Study of the effect of number of stringers ; Egle and Soder (19) 

studied the variation of the minimum frequency of a stringer-stiffened, 

circular cylinder with the number of stringers, keeping the total cross- 

sectional area (LA ) and the total torsional stiffness (L G J ) of the 
s s 

stringers constant. This is a reasonable approach for studying the 
explicit effect of the number of stringers. However, the implementation 
of "total" stringer properties being constant while the number of 
stringers is varied is more difficult in the experimental study than in 
the analytical study. 

In order to avoid this difficulty in the present study, the cross- 
sectional properties of all the stringers were assumed to be the same 
while their total number varied. Table XVII presents the variation of 
the natural frequencies of various circumferential modes of an internal 
stringer-stiffened freely supported elliptical cylinder with the number 
of equally spaced stringers. The geometric and material properties of 
the stringers are given in the footnotes of Table XVII. In order to 
visualize the variation of the frequencies of various circumferential 
modes with the number of stringers, some of the results of Table XVII 
are plotted in Figure 5. As is evident from Figure 5, the overall 
effect of the stringers is a lowering of the frequencies. This effect 
is greater on the frequencies pertaining to lower circumferential wave 
numbers. The rate of decay of frequencies due to the presence of 
stringers is greater for small numbers of stringers and diminishes with 
an increase in the number of stringers. 
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TABLE XVII 


STUDY OF THE EFFECT OF NUMBER OF STRINGERS ON 
THE FREQUENCIES OF A FREELY SUPPORTED 
ELLIPTICAL CYLINDER 3 (Hz.) 




e = 

= 0.760, m = 1 



\l c 

0 

- 2 

4 

8 

16 

n \ 






1 

1238.0 

1159.0 

1090.0 

984.5 

831.3 

3 

491.1 

470.4 

448.3 

450.2 

433.7 

7 

139.5 

121.1 

121.1 

122.7 

114.5 


183.5 

184.2 

184.3 

145.7 

141.9 

9 

226.6 

214.8 

212.7 

208.5 

204.9 

11 

263.7 

262.1 

256.5 

258.1 

224.6 

13 

380.0 

373.7 

368.5 

347.8 

290.6 

a) The geometry and material 

properties of 

the unstiffened 

shell are 

given 

in Reference 

(16). 





b) Circumferential mode number. 

c) Number of equidistant internal stringers. The properties of the 
stringers are: 


A s l 

- 0.1037 sq. in. 

z ls l = -°- 0475 in 

yy sSb 

= 0.005957 inf 

z... =.-0.2340 in 
2sl 

*zzsj l 

= 0.001285 inf 

y ls l = °-° in * 

^yzs l 

= 0 

y 2sA = °*° in * 


=912.5 lb. -inf 

00 

00 

m 

CM 

o 

o 

o 

• 

o 

II 

to 

Q- 


E . = 10.6 X 
si 

10 7 lbs .-sec f/ inf 
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0.01 X FREQUENCY (Hz) 



NUMBER OF STRINGERS 


Figure 5. Study of the Effect of Number of Stringers 
on the Natural Frequencies of a Freely 
Supported Elliptical Cylinder with 
s = 0.760, m = 1. 
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Ring- and stringer-stiffened elliptical cylinders : This section 
presents results for a stiffened, noncircular freely supported cylinder 
with large numbers of rings and stringers. The frequencies of the 
unstiffened freely supported elliptical cylinder with e = 0.760 are 
presented in Table XVHI. To study the effect of large numbers of ring 
and stringer stiffeners, 16 internal stringers and 11 internal rings 
were added to the above elliptical shell. The geometric and material 
properties of the rings and stringers are assumed to be the same and are 
listed in the footnotes of Table XVII. The frequencies and the mode 
shapes of this shell were obtained using the present analysis. Table 
XVIII presents some of the frequencies. Figure 7 shows some of the 
axial mode shapes and Figure 8 shows some of the circumferential mode 
shapes. To visualize clearly the effect of the large number of rings 
and stringers on the natural frequencies, some of the frequencies pre- 
sented in Tables XII, XIII, and XVIII are plotted in Figure 6. The 
results presented in Table XVIII were obtained with 5 axial mode func- 
tions and 6 even and 6 odd circumferential mode functions. It is quite 
evident from Figure 6 that the frequency curves of the ring- and 
stringer-stiffened shell under consideration, are more or less similar 
to those of the unstiffened shell; however, they are bodily shifted to 
the left. The minimum frequency of the stiffened shell is more than 
three times the minimum frequency of the unstiffened shell. The fre- 
quencies of the stiffened shell are consistently higher than those of 
the unstiffened shell. It should be noted that even though the ratio of 
number of rings to number of stringers in this problem is about 3:4, 
the effect of rings is predominant. Figure 6 reveals that the frequency 
curves for various m values tend to merge as n increases. The 


55 



TABLE XVIII 


FREQUENCIES OF 16 STRINGER 3 AND 11 RING 3 INTERNALLY 
STIFFENED FREELY SUPPORTED ELLIPTICAL CYLINDER b 
WITH e = 0.760 (Hz.) 


c 

n 

d 

m 


i 

3 

i 

741.0 

1703.0 

2 

444.9 

1303.0 

3 

437.9 

974.3 

4 

743.7 

973.5 

5 

1155.0 

1340.0 

6 

1868.0 

1998.0 

7 

2924.0 

2959.0 


a) The stringers and the rings have 
identical material and geometric 
properties which are given in the 
footnotes of Table XVII. 

b) The geometric and material pro- 
perties of the shell are given in 
Reference (16). 

c) Circumferential mode number. 

d) Axial mode number. 
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0.01 X FREQUENCY (Hz) 



CIRCUMFERENTIAL MODE NUMBER 

Figure 6. Comparison of Frequencies of Unstiffened, and Ring- 
and Stringer-Stiffened Freely Supported Elliptical 
Cylinder with e = 0.760. 
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a 



= 5, n = 5 ; 1739 Hz 


Modes 
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CIRCUMFERENTIAL 


Figure 8. 








CIRCUMFERENTIAL MODE m=l,-n=7; 2924Hz 

Figure 8. (Continued) 
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CIRCUMFERENTIAL MODE m=3, n = 4; 973.5Hz 



CIRCUMFERENTIAL MODE m=3,n = 5; 1340Hz 



CIRCUMFERENTIAL MODE m = 3, n= 6 ; 1998 Hz 


Figure 8. (Continued) 


6 . 



CIRCUMFERENTIAL MODE m = 3, n = 7; 2959Hz 


Figure 8. (Continued) 
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Conclusions 


1) There is weak circumferential modal coupling due to the presence of 
stringers in both circular and noncircular cylinders. 

2) The stringers contribute more to the total kinetic energy of the 
structure than to the strain energy. Therefore, the stringers tend to 
reduce the natural frequencies. 

3) The rings contribute more to the atrain energy than to the kinetic 
energy of the structure. Therefore, the rings tend to increase 

the natural frequencies. The influence due to the presence of rings 
is more than the stringers. 

4) Reasonably accurate tesults for ring- and stringer-stiffened shells 
may be obtained by considering the same number of circumferential mode 
components as are necessary when the stringers are not present. 

5) The reduction-of-frequencies effect due to the presence of stringers 
is greater on the frequencies associated with the lower circumferential 
wave numbers. 

6) The rate of decay of frequencies due to the presence of stringers 
is greater for small numbers of stringers and diminishes with the 
increase of number of stringers. 
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APPENDIX A 


DERIVATION OF THE COMPATIBILITY RELATIONS 


The compatibility relations of the stiffeners were derived based 
on the assumption that the stiffeners are attached to the shell along a 
line of attachment of infinitesimal width. This assumption is probably 
valid when the stiffeners are closely riveted with a single row of 
rivets. 

The displacement vector of any point in the cross-section of the 
i^k stiffener can be written (in vector algebra notation) as 

W - {'Vi } + (4 * fo/scil ■ 1 - { : £ -SUr. (Al) 

where q^ = The displacement vector of an arbitrary point in the 

cross-section of the stiffener; 
q = The displacement vector of the shear center of the 

stiffener; 

oj = The angle of rotation vector of the stiffener; 

R . . = The position vector of the point with reference to the 

i/sci 

shear center. 

These vectors may be expanded as follows: 


u 


sex 


SCI 


= < v 


sex I 




sex 


x = r,s 
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where (see for example, Reference (30)) 


W . a V . 
SCL.W SCX 

U xi R . “ R . 

sex sex 


x = r,s 


* " W - 

ox sex.x 


for rings 


v for 'stringers 

scs,x 


Also, (see Figures 2 and 3) 


R / 
r/scr 


R / 
s/scs 


where the vector components x' , y ‘ , and z are referenced to the shear 
center (sc). 

Substituting the above equations into equation (Al) , the compati- 
bility relations of rings and stringers result. 

For the rings: 


scr ,x 


s- W-kJ* 


scr ,x 
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For the stringers: 


/ / 

-z w - y v 
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(A3) 


scs 


scs 


'(■ 


, W A v 

SCS.O SCS 

y '-R— ^ - 5 — 


SCS scs 

Another set of compatibility relations were obtained to relate the 

shear center displacements of the stiffeners to those of the shell at 

the line of attachment by replacing r by scr, q r by q, z 7 by z ^ r > x 7 

by x , and R by R in equation (A2) and s by scs, q by q, z' by 
].r scr scs 

•z. > y' by y , and R by R in equation (A3). 

Is J. s scs 


For the rings: 


-z, w, 
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r \ _ f 
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j~ X lr 1 
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CD 
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l) + 

1 R U, 0 Z lr ' 

\ R R/ 


(A4) 

o ir \ k. R/ ^ 

x. w, 
lr ’x 

The cross-section of the ring was assumed to be symmetric with respect 
to the normal to the shell surface. Hence, the above equation reduces 
to 

' -z. w, 
lr x 


kJ-M* W£-i) 


(A3) 
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APPENDIX B 


ENERGY EXPRESSIONS OF RINGS AND STRINGERS . 


Ring energy functions: 
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)*■*.> 


d0 


scr x=x k 

(B4) 


Stringer Energy Functions: 


L E * 


U s (u.v.w) -£ -f !<*,!» ? x + + W ? => 

ext 1=1 0 * 9=9 . 


dx 


+ ~ - [<T-2 A .z .u, w, + { 
2 J > • s£ IsA x xx V. 


I *z, /, I ,Z. 

ZZS# lsjj + 2 zzsA lsl t „3 

R 2 R 


i}v? 

J XX 


{ 


I „z 2 I *2 
zzsA lsj& . ZZSjfc lsX 


a 

} V ’xx W ’: 


2 + -) , ’xx"’xx0 ' **Sj£“lS,fc"’xX 

L K 


I A 3 2 

+ A -z, .w. 


I ZZSJ0 Z 1SJ0 2 


+ R 2 W, xx0>> 


dx 


e=e 


E « p 

T*J< 


2 A .z_ u, w, +2 A 'z .z_ w. 


sjZ. 2s£ ’x ’xx sjJ lsX. 2s£ ’xx^>- dx 

0 9=9/i 




X XX 


-2 A . y !^ U ^ W, xx9 + ( A a + ^. ejls t 
si 2sl R \ Sj& 7 1s£ _ 2 


sA 7 ls£ r 2 


2 v 2 +2 f A y z . 

R / ’xx v sjrlsi, 1sj& R 


I .y, 

ime o' 1 


A sjj z lsl z 2sA y lsA 


+ A„ fl z_ y 


J v, w, +2 (- 
/ xx xx \ 


si 2sjj ls£/ ’xx ’xx 


I y 
yysjfc lsl 
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, A sl Z 2 s„ y ls^ v w Vs^lsA w3 + (tm sJl 

R / ’xx > xx0 2 ’xx0 ' R 

R 

K w, w, +w, w, )j> dx 
xx XX0 xx0 xx 1/^ 


i .y 


lsj& 


+ A „z, „z 


ls£ 


sj l lsi 2s A R 


0 = 0 , 




+2 


A sje 2 lsA y 2sA 

R 


U, x W, xx0 + ( 2 A s^ y lsJl y 2sji +2 


A sA z isjfc y i&e y 2&el 2 


} 


A sA Z lsje y lsA y 2sX, „ .. ^ A 


v > + 2 1 


’xx ’xx0 


s£ lsA-'2sX 


k sl Z lsi>2sl 

R 


> , w, 

xx 


A „z, „y 


sjJ ls£ J 2sA 


XX 


( w, w » ~ + w, ) /; 

xx xx0 xx0 xx// d 


dx 
e=e * 


+ ^.J<^-2 ( I y -g- s ^ 1 M + W z ^ y is^ + 2 

R 


0 

I _„z. 


yzsl 


+ i^lsi) +2 (izsspsl +2 


i y. 


- Vs^isd / W +w w )_2 IzsMllMiM w 3 \ dx 

R \ ’xx ’xx0 *XX0 ’xx/ 3 ’xx0 ^ 

R 0=0 


R' 


1 ™--> z isje y isje 


XX XX0 


l 

(B5) 


where: 


^zzsA ^zzcsjfc + A s£ y 2sX 


L . = I . + A .z,? „ 
yysj& yycsA si 2s l 


^yzs l *yzcs l + A sj& y 2si, Z 2sj& 


T (i.i.i) - itp , f<A y + I ,*? + (a , + iaai + »> 

s v ’ ’ ' *LSslJr^sl zzs£ ’x \ si „3 „ 2 / 

/It A K K 


1=1 0 
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i :. + i_,) ™, Q + 


— ( 

R 2 V 

+ (hzsL + !jssL) q+ a > aV > dx 

' *3 p 2 / *0 Si, ^ 0=Q 


r 2 \ yysj l zzs X' * 9 , yysX *x 

I * I 
\ 

R “ R' 


+ + 2 ^2 _ 2 ^ 
R 3 R 


a 

+ "si I <- 2 A si z l S / 5 -x + +2 


+ / I ZZSjj Z lsA +2 I zzs£ z is^ \ ^2 
’X V _2 R ' ’ 


A „z. 


R‘ 

I „z 2 . I .z. . 
zzsX lsX . zzsX lsX 


x 


+ 


R 


) o 0 

V, x W ’x 0 


2 + k si, Z lsJL 

' _ p. 


R 


\ 0 o 

) vw, 


2 o2. ZZSX lsX ..2 

- + A .Z, .W, + W, - 

G si 1 st ’x r 2 X 0 


A ,z 


sX lsX o 2 \ J__ 

+ w, e p> dx 


R 


0 = 0 , 


+ p , f<<2 A ,z, ,41. +2 (^M + A .«‘lM*2.i ) v 

H sX J sX 2sX ’x V R 2 / 

IS- 


A ,z n . A .z n .z„ . v A „z, „z„ „ 

2sX „ sX lsX 2sX ) oe. „ s X lsX 2sX » 2 

-rZ J VW, fl -rZ W, A 

r2 / e R 2 0 


/A „z 

- 2 (^ir 


+2 A sX Z 1 sX Z 2 sX W, 3 x/ > 


dx 


e=e, 


j < 2 (- 


A sX Z 2 sX y lsX' 


A sX y lsX + R ) UV ’x ' 2 UW 


A sX Z 2 sX y lsX oo 


R 


*x 0 


+ ( A v 3 Vsxyisx A sX Z 2sX y lsX y. 2 ,j£lsi. 
X BX y isX R 2 r ' ’X r3 

+l(k z V I yyso y lsX A sX Z lsX Z 2 sX y lsX 

V sX lsX y lsX R R 
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+ k sJl Z 2sZ y lslL/ v ’x"’x 


0 . ( yysl y lsl . si 2sl y lsl\ „ . 

J v, w, +2 ( " + ) v, w, 

/XX \ 2 R / X 


x0 
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R' 


R 


R 


+ 


A sl Z lsA Z 2sl, y lsl, 

R 


Woo . o o W A si, y .. 
M W ’x W ’x0 +W ’x0 W W + _ W > 


S lsjft 0 2 
2 ' 0 


SJ& y Ls£ ( a o ,6o j 

+ — — V w > e w + ™» e A> dx 


0=0 


l 


+ p 


si 
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A sA y 2sl + 


A - Z ls/ y 
R 


sl^lsl^s! \ o. 

) uv s 


+2 ■t s&sfZls k ^ 

x R x9 


+ [ 2 A ,y, ,y, . +2 , 2 


0 


si ls2 J 2s! 


R 


R 
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t» 2 / X „ 2 ' Sx» 1SX» 


K,y-,, 0 y-. 


A sl z lsl y 2s2 

R 
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R 
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si lsl 2sl 


+ 1 , + txsspsi)*, i, +2 (. 

yzsA R /xx' R 


yzsl 
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sl y 2sl 00 
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-2 


A sl y lsl y 2sl 00 Y A ,l Z Ll y 28l 
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R 

I .z 


R 


yzsl lsl \ / o o ,o o N „ yzsi y isl lsl » 2 

+ — 1 — ^ • J V w, w, . + w, -w, J -2 — 1 1 w, . 

R /' x x8 x0 x/ „ 2 x0 


A sl y lsl y 2ss, , A s! y 2s! , . » W 

+2 — — L w?q .+ — — ^ww„ e + w, Q *)^> dx 

K R /n *\ 


R 
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APPENDIX C 


MATRIX ELEMENTS AND INTEGRALS 


The matrix elements of Equation (32) and the circumferential and 
longitudinal integrals involved in these elements are presented in this 
appendix. The closed-form expressions for the longitudinal integrals 
were obtained with the help of a table of formulas for integrals derived 
by Felgar (31) . The circumferential integrals were evaluated numeri- 
cally using the 8 -point Gaussian quadrature method with four subintervals. 

The elements of the mass and stiffness matrices of a ring- and 
stringer-stiffened noncircular shell may be written as follows: 


Contribution of the Noncircular Shell 

A -- = S, IS1 IX a + (S 3 IS1 3 + S 3 IS1,) nn IX 3 
mn,mn l a J 3/ 

D — — = S 4 nISl 5 IX 3 — S 2 nIS 1 = 1 X 2 

mn,mn J “ 0 

E mn,mn = S 4 IS 1 s IX 3 -S B ISl s IX i + S 3 nnISl 7 IX 3 

™ = SxnnlSlelXg + (S 3 IS1 9 + S 0 IS1 3 )IX 2 + S 5 IS2 IX 5 

mu y mil 1 

F mn,mn = S 1 nIS1 B H 5 " S,nISl B IX 4 + S^IS1 3 IX 3 - S 7 IS2 3 IX 4 
-S 5 IS2 3 (l-n 3 )IX s + S^ISZilXg 

C mn,mn = < S i IS1 e + S 5 IS1 4 )IX 5 + S^ISl^ + (n 3 n 3 - n 3 -n 3 )ISl ,IX 5 ! 

-S 7 ISlg (ii^Xg + n 2 IX 4 ) + S B niilSl 3 IX 3 + S 5 jnnIS2 i + (nil 3 
-n) IS2 3 + (n. 3 n - n)IS2 4 J IX 5 - S, (nIS2 5 IX 3 + nIS2 3 IX 4 ) 

N mn,^= 2 

V,mn = 2 P c hISl 9 IX s 
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S -- = 2 p hisij IX 5 
mn,mn c 


(Cl) 


where IS^ to IS2 3 are circumferential integrals, IXx to IX 5 are 
longitudinal integrals, and S x to S 8 are constants defined in Appendix D. 
The circumferential integrals are defined as follows: 

TT 

ISli = Jr cos n9 cos n9 d9 
0 
TT 


IS1 2 = J ^ sin n0 sin n© d9 


0 

IT 


IS1. 


r i 

= I sin n9 sin n© d9 

J r 3 


0 

TT 


IS1 


i R 3 


cos n9 cos n8 d0 


0 
TT 

ISl s = J cos n9 cos n9 d6 
0 
TT 


IS1, 


IS1. 


IS1, 


= J si 


in n9 sin n9 d9 


0 

TT 


= I ~ sin n9 sin n9 d9 
J R 2 


0 

TT 


-f I 

J R 


cos n9 cos n9 d9 


0 

tt 


IS1 9 = J R sin n9 sin n9 d9 
0 
TT 

IS2 X = J ^ {( r ) } s ^ n n ® sin 


d0 
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IS2 3 



sin n0 cos n0 d0 


TT 


IS2 3 



sin n0' cos n0 d0 


IS2 4 


fl(I) 

J „ 2 NR/, a cos n0 sin n0 d0 

a R 0 


TT 

IS2 5 = J GL cos n0 sin n0 d0 
0 ’ e 


(C2) 


The matrix elements of the antisymmetric mode equations for the 

shell are identical in form to the above equations and are obtained by 

interchanging Sine terms with Cosine terms and vice versa. Furthermore, 

(m) must be replaced by -(rO . It was found that if the cross- 

XR/ >e 'R-se 

section of the shell is symmetric with respect to the horizontal axis of 
the shell, there is no coupling between the even and odd terms of n and 
n. Thus, in the analysis of elliptical cylinders, two computations must 
be made in both the cases of symmetric and antisymmetric modes (with 
respect to the vertical axis); one with all even terms of n and n, and 
the other with all odd terms of n and n. 

The longitudinal integrals may be defined by a general axial mode 
f unc tion 


as follows: 
IX 

l 


IX 2 



§- dx 
m 


0 



§ ' § - dx 
m m 


0 


§ 

m 
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IX, 


IX, 


IX, 


-U 


0 

a 


-j 


0 

a 


-i 


$- dx 
m m 


$ §- dx 

m m 


§ f - dx 
m m 


(C3) 


Substituting Equations (28b to 28e) into the above equations, the 
longitudinal integrals for various boundary conditions may be written 
as : 

For freely supported cylinders: 


IXi 

m rr 

1 

l 



2a 3 

3 2 ( 

1 


IX S 

= -IX 3 

II 

It 

ro 3 
61 =. 

1 For 

m = m 

IX 5 

a 

1 



= 2 

l 



IX 1 

to IX 5 

" ° • 1 

i 

| For 

m ^ m 


For clamped-free cylinders: 

IX ! = 


p 4 a 

m 


IX a = 


Of P (2 + a P a) 
mm mm 


^ to#’ - af |) 


m m 


-P-8 (a 8 - a-8-)l 

mm mm m m J 


aj.( 2 - a 8_ a) 

mm m. m 


m = m 


/ m 


m = m 


m ^ in 


m = ffl 


IX, 


48! (a-P- - o: P ) 


m m m m m 


m m 


rf ( - 1)m+ " p « +e a 3 ] ”’ ,a 


(C4a) 
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a _P (2 - a 8 a) 
mm mm 


m = m 


ix 4 - < 48^(0 0 - cy_ 3 _) 

J m mm mm 


m mm m m' 

3 4 - 

m m 


[<- 1 , 


m + m 


P- 2 + P^l 

m m J 


IX, 


a 

0 


For clamped-clamped cylinders: 


3> 

m 


* IX 


IX 2 = -IXo = -IX 4 = 


a 3 (aj* a ” 2 ) 

m m m m 

48ffe a (a_P_ -a 3 ) 

mm mm mm 

3 4 - 3! 

m m 


m ^ m 

m = m 
m ^ fh 

m = m 

m ^ m 
m = m 


[<-i> m+i + i] 


IX, 


.0 


m = m 


m £ m 


For free-free cylinders: 
m = 0 

IX-l = IX 3 = IX 3 = IX 4 =0 


IX 5 = a 


IX i = IX 3 = IX 3 = IX 4 = IX 5 = 0 
IX x = IX 3 = IX 3 = IX 5 = 0 


IX 4 = 


4or- ,3- , 

m-1 m-1 


m = 1 


IX X 

II 

0) 

X 

M 

II 

II 

CO 

X 

w 

s 

11 

K 

u> 

II 

0 

IX x 

= IX 3 = 

II 

X 

M 

d 

IX 3 

1 

a 5 

ix 5 = 

a 

12 

IX ! 

11 

M 

X 

u 

II 

1 X 5 = 

0 


m = 0 

m = 1 
m S 2 

m s 2 even only 
m > 2 odd only 

in = 0 
m = 1 

m > 2 


(C4b) 


m / m 


(C4c) 
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IX 2 = 


IX 4 = 


4 

a 


VO 


- - 2a- .P- , 
a m- 1 m- 1 


V 0 


2 2 


ix i = IX 2 = IX 4 = IX 6 = 0 


4p .a , 

m-1 m-1 


IX, = 


IX, = IX 4 = IX = = 0 


IX 2 = 


IX, = 


4 

a 


- - 2a ft 
a m-1 m-1 


IX = 
1 


m > 2 odd only 
m s 2 even only 
m > 2 odd only 
m 2 2 even only 


m 2 2 even only 
m > 2 odd only 

m > 2 odd only 
m s 2 even only 
m > 2 odd only 
m 2 2 even only 


IX a = 


'P 4 ,a 

I m-1 


or iP„ i( a m if* i a + 6 ) 
m-1 m-1 m-i m-1 


4 MlMlM-lMl" o; m-lMl ) 

Ml - Ml 


m jt m 




D 


+ (-D 


m + m - 2 


] 


O' l^m l^ 2 " 1 

m-1 m-1 m-1 m-1 




m ^ m 


m = m 


IX, = 


Mi - Mi 


l^m-fm-1 m- 1 m - V r 
q4 q 4 L 

1 


m / m 


+ (-D 


m + m - 2 ’ 
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a m 1^ 2 ~°L I 3 ) 

m -1 m- l m -1 m -1 


m = m 


IX, = 


IX 5 = 




1 m- 1 H m - 1 m- l p m- 1 ' L 

Cl - Cl 1 ' 

+ ( _l )m + m- 2 j 


m ^ m \ m s 2 


m = m 


m / ffl 


(C4d) 

The longitudinal integrals in Equations (C4a, C4c, and C4d) vanish if 
m+m is odd and are nonzero if m + m is even. 


Contributions of Stringers 


mn.mn 


-I K 


1=1 

L 


cos n 0 cos n9) 


T’ /! oa 2l v _ v 

D __ = )\( -SS + —5 SS 30 — ) cos n0 sin xi© \ IX : 

mn,mn n K K / s q~q 


0=0 


i 

SS 2l \ 


1=1 
L 

E __ = ) / - SS , cos n0 cos n0 -SS 9 cos n0 cos n0 

mn,mn 

Z= 1 

SS : 


, _ cos n 0 sin nS 21 _ » . ^ , v 

+ SS l3 n — - — n cos n9 sxn n0 X IXi 


/0=E 


b .. -v<(ss 3+ ^ 

mn,mn 4 j\\ p s 


1=1 


SS 5 ss 6 

o + ~ 

R 


SSi 4 SS 1 5 


+ nr + SSi3 + ~ r + SSs2 


R 


C ] 


SS 23 SS 2s SS^ \ 

+ — - — - — - — - ) sin n0 sin n0 ^ IX x 

R R r3 / 


+ (GJ) , Bl ° " 6 Sl ° 56 IX 3 | 

Sl . "R 3 0=0 


Z 
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inn, mri 


£=1 


^ 87 ^ _ / o:5 i 

= ) \ ( 1 — — ) n sin n9 sin n0 . + (SS ls - — — 

mn \> 2 R/ ' R 

K 


ss,, ss 


1? aa 18 


R 


\ /SS]^ SS g v 

+ SS J sin n0 cos n0 + - — — — J ft sin n9 sin ii0 


R 


SS 


24 


R 

SS; 


R 


n sin n9 sin n.0 + ( SS 21 + 


/ sb 2S\ 

\ SS 21 H R J si 


in n9 cos nS 


/ SS 2 gv 

+ ^ SS 2 s H J si 


in n9 cos n© 


✓ SSqq SSg^v 

4 — - 1 In s in n© s in n9^^ IX x 

' R R a / 

<H e ^ sinn :r fi9 > c 


9=0, 


'0=0, 


Lt 

l ,inn = 


SS, 


Cmn.mn = )< SS 3 cos n0 cos n0 + SS e cos n9 cos n0 H nn sinn0sinn0 

2 

4= r 


R‘ 


SS 14 /SS^ r, 

• + SS 1Q cos n9 cos n0 + nn sin n9 sin n0 - ( — - — 

_ 2 ' R 


SS 


18 


R 

SS 25 

R 

SS 3 g 


R 


)( fi 

(* 

(“ 


cos n0 sin n9,+ n sin n0 cos n9. 


cos n9 sin n0 + n sin n9 cos n0 


cos n0 sin n9 + n sin n9 cos n0 


SS 


37 - - \ 

nn s in n9 s in n© > IX 3 

R 3 ^0=0^ 


<« 


. , - sinn9sinn0\ 

+<(CJ)sinn __ y a . 
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N 


mn 


- -i<*. 

,nm 4 j\ A 


1 = 1 
L 


NN 


ran 


.•■K( 


cos n 0 cos n©X IX 2 

X e = 0 je 

Tl7 ^-30 \ 

-Tj s + — T 29 + j cos n6 sin n 0 


4=1 

L 


)> IK -- 
/ 0=0, 


vs -K- ( t * + Tia ) ” s » e ■=» s 5e + (- 

T 3 °\- - V. 

- — — )n cos n 0 sinn 8 J> IX 2 

R ' A A 


17 

R~ 


4=1 


r-> // T 7 T 8 T 19 

,mn =L\S T * + - 2 + -T +T '* + -^ ' 
4=1 R R 


V 


9= V 


Tl 9 T 


20 


l 32 i 33 


— + t 31 + -^ r 


l 34 


R' 


\ - X y7 T 4 Tg T 

) sinn 8 sin n 9 > IX 2 4 -\( Ti H -i f — 

/ / e = 9 . X r3 r2 R 


e 




R R 2 R 2 R 2 


T13 T 14 T 18 T 31 \ _ 

H — — 1 H 1 J sin n8 sin n0^^. 


IX, 


0 = 0 , 


R 


mn 


,« =I<( 


4=1 


T rp ¥¥i iti fn rii 

5 *-12 *-S 1 1 3 J -15 *-25 

R 3 R 2 R R r 2 r 2 


*^ 39 \ _ 

+ - — ) n sinn© sin fi0 

_ C> / 


in fi0 \ 

^ 0=0 


Tn T 8 Ta 4 

2 


// T l 1 T 8 

1X8 + \(tt + t + 


R ' 


*-2 0 

l 32 

*-33 

* 38 X 

R H 

h R ‘ 

R 



:in n0 

/ o=P 


sin n0 sin n0 X IXa 

0 = 0 , 


4 


<(' 


T 32 t 


23 


L 3 5 


+\l T 2i R g— + T 17 + T 3 o + -g— + T 30 


-3 7 \ 

T” / 


sin n0 cos fi0 


'**1 V 

T 38 \ V 

+ — J sin n 0 cos n(T> 0=0 IX 5 
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s 

mn,mn 


L 

= £<(t 3 + Tg + T 14 ) cos nG cos nG^ IX 2 +<^T 4 + T s 

i= 1 ' 9=9 * 

+ T x 4 + T ls + T 31 ) — sin nG sin nG IX 5 

' P 3 / Q = Pl 


+ ^Tx cos nG cos nG 


S ix 5 +<Tt 7 + t iS 

/ e=e. 


T 

T 34 )— sinnfl sinn6'\ IX 2 

' r* 3 n _rs ' ' ** 


£)<= 

<0 


0=G, 


cos n0 sin nG + n sin nG cos nG ) IX S 


1 


+ \ V T 40 + T 41 ~ T 3 S " T 2 7 


)(« 


cos nG sin nG 


+ n 


sin nG cos nG 


)> 


IX; 


0 = 0 , 


(C5) 

where SS^ to SS 30 and Tx to T 42 are constants defined in Appendix D. 

Contributions of Rings 

K 

A — = ) C xn 2 fi 3 IRlx X , + C xnfilR4i X x + C x(IR4 s nn 2 +n 2 ii IR4 3 )X x 
mn,mn U A A A A ox 

k=l 

4- C -2 1 nn IRI 5 X 1 
K 

E mn fin = Z C m^RlgX x- C 4 n 2 fi 2 IRlx X x+ C 1 nIR4 4 X x - C 4 nn IR4xXx 
k=l 

; * C 4 n n IR4 3 X x - C 4 nn 2 IR4 2 Xx d - C 2 x*^ IRlgXx 
- C 3 5 nn IR1 5 X x 
K 

B -- = ) ICetin IR1 3 +C 3 nn IR 1 x)x 2 + (c 5 IR2 X + c 0 IR2 2 + C 7 IR^ 
mn,mn o v / \ 

k=l 
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+C e IR2 4 )nflx 3 -f(c u IRl s + C X 2 IR3 X + C X 3 IR3 2 ) nnXs 
+ C .3 IRA 1 Xs + C 3 (lR4 2 n + IR4 3 n) X 2 + ( C 5 IR5 18 


+ C a {lR5 x + IR5 2 . } 4- C 8 {: 


IR5 3 4- IR5 4- C i g IR5 X 2 
4 - |c 9 ( nIR5 7 + nIR5 s )+ C 19 (nIR5 9 4- nIR5 e )+ C 5 (nIR5 X0 
+ nIR5 1 1) + C e (nIR5 X 2 4- nIR5 13 ) 4- C 6 (nIR5 xs + nIR5 X4 ) 

+ Cg^nIR5 X 7 4- nlRSj.g^j’Xjs 4- £c x2 IR 6 X 4-C X7 IR 6 2 
+ C 13 IR 6 3 ]-X 2 4- jc 17 (klR 6 4 4- nIR6 5 4-nIR6 xo + nIR 6 xx ) 
4-C x 2 fnIR 6 e 4- nIR 6 .^4- C x 4 ^nIR4 4 + nIR4s ^4* C 2 q 


nIR 6 f 


-f nIR 6 9 ) f Xg 


)} 


K 


k=l 


F mn,f5n = Z<^ C 3 nfi3 IR1 i + C s nIRl 3 4-(c 8 IR2 4 + C 9 IR2 3 4- C 5 IR2 X 

+ C 6 IR2 2 ) nn 2 +• C 9 IR2 3 n + (c X 4 IR1 2 + C x 5 IR3 2 
+ C 12 IR3 1 )nn a +(c 14 IR1 2 4- C 16 IR3 2 )n 4- C 3 n 3 IR4 2 
4" C 3 IR4 3 nn 4- C 3 nXR4 x 4- C 8 XR3 } q 4* C XR5 x 2 4-XR5 x 4 ^ 

4-C e IR5 x 6 4- C x 9 IR5 8 4- {c 9 IR5 7 4- C 5 IR5 XX 4- C X 9lR5 9 
4 - C 6 (lR5 13 4- IR5 X5 ) 4- C S IR5 X 7 } nn 4-{c 5 IR5 xs 4- C s (lR5 x 
4- XR5 2 ) 4- C 8 (lR5 3 + IR5 4 ) 4- C X 8 IR5 5 ]- n 4- C 9 IR5 6 
4- {c X 2 IR 6 6 4- C 2 o IR 6 8 4-C xa (lR6 xo 4- IR 6 4 )} n s 4-£ 17 (lR6 5 
4” IR 6 x x ) 4- C X 2 IR 6 7 4" C X 4 IR4 8 4- QsoIR^g^ nn 4 -^C X 2 lR 6 x 
4- C X 7 IR 6 3 4- C 12 IR 6 3 ]- 5 4- C x s (lR 6 4 4- IR 6 xo ) 

4- C x 4IR4 4 ^^ X : 2 
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K 
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where IRli to IR 6 ]! are circumferential integraXs and Xi= § , $I| 

A A A m m 1 x=x^_ 

and X ? = $ $- and C, to C 40 are constants defined in Appendix D. 
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The quantities X x and X a for different boundary conditions are 

defined as follows: 

For freely supported cylinders: 

- 3 mTT x, mn x, 

n mmTT k k 

X = 2 cos cos 

* ^ a a 


mn X, mn x* 

_ A t K 4 

X 3 = 2 sin sin 


(C8a) 


For clamped-free cylinders: 

Xl = P m P m { Sinh P m*k + Sin P mV “m ( COSh P mV COS P m X k)}{ Sinh P S*k , 
+ sin P-^ - Qf- (cosh P cos P-xJ} 

X 2 = {cosh P^ - cos P^ - a m (sinh P^ - sin P m xJ}{cosh 
Pi*k • “m ( sinh P i*k ' Sin P m X k)} 
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For clamped-clamped cylinders: 

Expression is same as clamped-free but a ' s and 

3 's are different, 
m 

For free-free cylinders: 
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APPENDIX D 


CONSTANTS OF MATRIX ELEMENTS 


This appendix contains the constants used in equations (Cl, C5, 

and C6) of Appendix C. These are various combinations of the stiffner 

properties given in the list of symbols. 
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12D(1-V) 
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s 4 

S s 
S s 
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ss 2 

ss 3 
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ss 5 

ss 6 

ss 7 


D(l-V) 

24DV 

h 2 

2D 

3D(1-V) 

2D V 

4D(l-v) 

E sl A sl 

E sZ 1 zzsJl 

E ,1 
si yys l 

E sl A sl Z lsl 

E sl*zzsl Z lsl 

^ E sl^zzsl z lsl 

E sl^zzsl z lsl 
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ss e = 

SSg = 

55 10 = 

55 1 1 = 
SS 12 = 

SS 1 3 = 

5514 = 

5515 = 
SS 1 8 = 

SS 1 7 = 

SSi 8 = 
SS 1 g = 

55 20 = 

5521 = 

55 22 = 

55 23 = 

SSg 4 = 

SS25 = 
SSg a = 
SS 2 7 = 
SS 2 8 
SS2 9 


E s£ A s£ Z 1sj& 

E si A sZ Z 2si 
2 E si k si Z lsZ Z 2sZ 
E sZ A sZ y lsZ 
E sZ A sJL Z 2si y lsl 
E sZ A sZ y lsZ 

FT v ^ 

si yysl y lsl 
2 E sZ A sZ z 2sZ y lsZ 
E sZ A sZ y lsZ z lsZ 
E s£ I yysji y lsj£ 
E si A sZ Z 2sZ Z lsl y lsZ 
E sZ A sZ Z 2sl y lsl 
E sZ A sZ y 2sl 
E sZ A sZ y 2sZ z lsZ 
2 E sl A sZ y 2sZ y lsJL 
2 E sjJ A s£ y 2sA y ls£ Z lsA 
E s^ A s < 0 y 2sJi y is < C Z lsjl 

E s/s* y 2sW s * 

2 E sX I yzsA y lsjl 
2 E sZ^yzsZ y lsZ 7 "\.sZ 
E sZ 1 yzsZ 
E sZ^yzsZ Z lsZ 
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c 2 

= 2 

E.A, 
rk rk 

C3 

= 2 

rk xxrk 

C 4 

= 2 

E rk*zzrk Z lrk 

C s 

» 2 

E.A, zf . 
rk rk Irk 

C s 

= 2 

FT 2 ® 

rk xxrk Irk 

C7 

= 4 

E rk A rk z lrk 

C 8 

= 4 

E , I . z. , 
rk xxrk Irk 

Cg 

■ 2 

E , A z, , 
rk rk Irk 

C10 

= 2 

n rk zzrk Irk 

Cu 

= 4 

E rk A rk Z 2 rk 

C 12 

= 4 

E rk A rk Z lrk Z 2 rk 

C 13 

* 8 

E rk A rk Z lrk z 2 rk 

C 14 

- 2 

E rk A rk Z 2 rk 

C 15 

= 6 

E rk A rk Z lrk Z 2 rk 

Ci a 

= 2 

E rk A rk Z lrk Z 2 rk 

Ci 7 

= 4 

E rk A rk Z lrk Z 2 rk 

C 1 8 

= 4 

E , I , z5 
rk xxrk Irk 

C 1 9 

= 2 

E ,1 , z, , 

rk xxrk Irk 

C20 

= 2 

E rk A rk Z lrk z 2 rk 

C21 

= 2 

< GJ >rk 

Cas 

= 2 

^rk A rk 

C33 

= 2 

^rk^zzrk 
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Cs 4 = 2 p , I , 

* K rk xxrk 

C 2S =2 CGJ) rk z lrk 
C 3S = 2 (GJ)^^ 

C 27 =4 (GJ) rk z lrk 
c 28 =4 P rk A rk z i rk 
C39 - 4 P rk I zzrk z i rk 

c 3° - 2 P rk A rk z i rk 

C31 =2 P rk A rk z lrk 

c 32 =2 

C33 - 4 

C 34 = 2 P rk I zzrk Z lrk 

C 3 s =4 P rk A rk z 2rk 

C 36 = 2 Prk A rk Z 2rk 

37 = 4 Prk A rk Z lrk Z 2rk 
C 3 8 8 Prk A rk Z lrk Z 2rk 

c 39 = 12 P rk A rk z i rk z 2rk 
C 4Q =4 P rk A rk z lrk z 2rk 
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